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ABSTRACT 


Forty years of precipitation and surface temperature data observed over 
261 Local Climatic Data (LCD) stations in the Continental United States was 
utilized in a ground hydrology model to yield soil moisture time series at each 
station. Using Thornthwaite's (1948) empirical functions to calculate potential 
evapotranspiration and Nappo's (1975) moisture availability function to then 
calculate actual evapotranspiration, a month- by-month soil moisture dataset 
was constructed for each year. The data was averaged to obtain forty year 
mean fields. The forty year monthly mean soil moisture field shows a marked 
seasonal cycle across most of the country, with driest conditions ocurring 
in September. The calculated soil moisture is close to the saturation 
limit throughout the winter months. 


1) Introduction 


Several papers have mentioned the role of soil moisture as an important 
boundary forcing parameter influencing the seasonal mean circulation and 
convective precipitation (eg. Sutcliffe, 1956, Namias , 1959 Schlckendanz , 1976). 
Some recent investigations with General circulation models have been carried 
out to study the influence of soil moisture on precipitation and climate 
(e.g. Walker and Rowntree (1977), Shukla and Mintz (1982), Miakoda et al . 
(1979)). However, to date there has been no systematic direct measurement 
of soil moisture as a meteorological variable. Therefore researchers have 
either created their own soil moisture fields via hydrological models 
(e.g. Mlntz and Seraflni, 1983), or inferred soil moisture values from satellite 
data (e.g. Soer, 1980, Schmugge et al . , 1980). The authors themselves faced 
this problem in the following study of the influence of soil moisture on 
ensuing precipitation. For this purpose, observed monthly mean surface 
temperature and precipitation data supplied by the National Climate Center 
(NCC) was used in a ground water balance model to calculate soil moisture 
time series at each of the 261 local climatic data (LCD) stations across 
the contiguous 48 states in the United States. The mean fields resulting 
from this calculation are approximately representative of climatic normals. 

Forty year monthly means and standard deviations of temperature, precipi- 
tation and soil moisture over 261 NCC LCD stations are presented on U.S. maps 
as well as in tabular form on microfiche. The raw data utilized is described 
in section 2, and the surface water balance model used to derive the soil 
moisture from the temperature and precipitation fields is described in 


section 3 


Monthly Local Climatic Data (LCD) consisting of precipitation and 
surface temperature fields, from 2.61 stations over the contiguous 48 states 
in the U«S» was obtained from the National Climatic Center (NCC) in Asheville, 
North Carolina. These monthly means were tabulated from actual data collected 
at hourly or three hour intervals. The data cover a period of 40 years 
(1940-1980). Incomplete data years at any given station were ignored as opposed 
to filling in the data by interpolation. Table 1 shows the station names, 
sequence numbers, WBAN identification numbers and station altitudes. The 
sequence numbers are shown in Fig. 1 to aid in locating the stations on 
the plots. Fig. 2 shows the number of complete data years (maximum of 
forty) at each station. 

Means and standard deviations of temperature and precipitation were 
calculated at each station for all complete data years. These fields are 
shown plotted on I’.S. maps in figure s*;:;» CH»u) and (4a, b) as well as in 
tabular form on the microfiche inside the back cover. 

3) Calculation 

a. Determination of evapotranspiration : 

In order to determine the evapotranspiration a ground water balance 
model similar to that employed by Mintz and Serafini (1982) was used. 

While this model is a simple representation of the extremely complicated 
processes of evapotranspiration on land covered by a variety of vegetation, 
it was considered adequate to reproduce the first order soil moisture from 
the temperature and precipitation fields. Since the model employs a moisture 
availability parameter based on observations over grass covered land, (Davies 


ORIGINAL PAGE IS 
OF POOR QUALITY 


and Allen, 1973) the calculation may be somewhat inadequate for dense forests* 
We expect It to be more accurate for drier regions with sparse vegetation* 

For a justification of the use of this 0-curve for all types of vegetation 
cover, see Sud and Fennessy (1982)* 

The potential evapotranspiration, Ep at each station was calculated by 
Thornthwaite's (1948) method, which was derived from observations of water 
use in Irrigated regions* The method yields Ep as a function of surface 
temperature, Ts, time of the year and latitude* The influence of different 
climatic regions on evapotranspiration is taken into account by introducing 
a yearly mean surface temperature parameter. The empirical relations are 
as follows.', 

(i) For Ts > 26.5°C (See Willmott, 1977), the following relation is used: 



Ep * a + b Ts + c Ts 

(1) 

where 

a - -415.8547 



b - 32.2441 



c - 0.4325 


(ii) 

For Ts < 26.5°C (Thornthwaite, 1948), 

the formulation is an follows: 


Ep - 16 (10Ts/H) k 

(2a) 

Here H 

12 

- 1 (Tsi/5) 1 * 514 

i-1 

(2b) 

and k ■ 

« a' + b'H + c'H 2 + d"H 3 . 


where 

a- ■ 0.49239 



b" ■ 17.92 x 10“ 3 
c" - -7.71 x 10-5 
d" - 6.75 x 10-7 


(iii) 


For Ts < 0°C 
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Ep = 0. (3) 

Equations (1) thru (3) yield potential evapotranspiration as a function of 
surface temperature for different surface temperature conditions* The heat 
index* H* accounts for differences in evapotranspiration in different climatic 
zones* Clearly it is a function of mean annual surface temperature* At 26*5°C, 
the potential evapotranspiration can be determined by either (1) or (2)* Obvi- 
ously the relations are incompatible at this boundary. However* this incompati- 
bility introduces only a small error. In an example equation (1) yields 
Ep - 4.50 mm/day and equation (2) yields Ep between 4.44 mm/day and 4.60 mm/day 
depending upon whether Ts was 26.5°C averaged over the year or 26. 5°C during 
the month in question and 0°C for the rest of the eleven months. 

Ep is corrected for the days in the month, D, and number of daylight hours, 
L, in the following way 

Ep - Ep ( D) ( L) (4) 

c 30 12 

The actual evapotranspiration, E is then obtained by multiplying Ep c by 
Nappo's (1975) moisture availability function, M as adopted by Mintz and 
Seraf ini (1983). The set of equations are: V 

E M Ep 

_ aw 

M « [ 1 - e w *] 

Here a is an experimentally determined constant. Its value was set at 6.80. 
w is the actual soil moisture in mm. The saturation field capacity, w*, was 
assumed to be identically equal to 150 mm over all land Irrespective of soil 
type or vegetation cover* 


(5) 

( 6 ) 






b. Determination of toll moisture: 




Neglecting all water transport within the soil, the soil moisture 
tendency equation is: 


iH ■ P - E 

at 


(7) 


where P represents the precipitation rate. 

Substituting from (S) and (6) we get 

-aw/w* 

in - P - l 1 - e J Ep (8) 

at c 

Obviously (11) is non-linear in w« 

It was solved by numerical integration using time steps of one day. 

In order to avoid precipitation input from the next month, the mean monthly pre- 
cipitation and surface temperature values were assumed to be constant for the 
entire month. Correspondingly Ep c was also calculated from (4) and held con- 
stant for each month. The soil moisture, w 1 *^ on day n+1 was obtained from its 
value w n on day n by a two step predictor-corrector technique: 


w n+l « w n + dw n At 
dt 


and then. 


(9) 


w n+l > w n + A 


dw n + dw 


irt-1 - 


dt 


dt 


At 


( 10 ) 


This calculation is fairly accurate for small time steps. The limiting con- 
straint, w 1 * 1 < w*, was Invoked after each time step. In other words if w n+ l 
exceeds w* at any time, the difference is used to produce runoff, the soil 
moisture time series was obtained by the above described calculation for 
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each station. For each year the integration was always started from the 
equilibrium soil moisture on January 1 obtained separately by cyclically inte- 
grating the forty year mean monthly data. This avoido complications arising 
from years with missing data and eliminates the build-up of any systematic 
errors in the calculations. 

The mean monthly soil moisture values and the corresponding standard 
deviations were calculated from the time series at each station. Plots of 
these fields are shown in Figure sets (Sa) and (5b) respectively. In addition, 
these values are tabulated by month and station on the microfiche appended 
at the end of this report. 

Discussion 

The 40 year mean monthly surface temperature fields are shown in Fig. 3a. 

The seasonal cycle is clearly depicted in these maps with July being the hottest 
month. The calculated potential evapot ranspi ration field (not shown) also 
reached its maximum in July. The standard deviations of surface temperature 
(Fig. 3b) were found to be larger inland and during the winter months, as would 
be expected. 

The monthly precipitation fields (Fig. 4a) also agreed well with the known 
climatology. The standard deviations of precipitation (Fig. 4b) were of the 
same order as the field itself. The monthly mean temperature and precipitation 
fields are presented in this report because the soil moisture time series was 
derived from their time series. The mean soil moisture fraction field shows 
a distinct seasonal cycle across the contiguous 48 states. 

During winter there are two regions with large soil moisture gradients. 

The mid and northern parts of the west coast are completely satuarced, as is 
is the majority of the country east of 95°W» However, the region between 115°W 
and 100°W is quite dry. The strongest gradients lie in the 100°W to 95 6 W region. 
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As spring progresses these gradients tend to spread out, with the entire 
country generally drying. By May very few stations across the country remain 
saturated. 

By late summer the soil moisture has generally reached its driest limit 
across the U.S. During September the entire country west of 105*W is extremely 
dry with soil moisture fractional values below 0.1. The rest of the country 
is also quite dry, with Florida being the only region with values above 0.5. 

As fall progresses into winter a gradual increase in wetness is noted. 

The eastern U.S. again becomes quite wet as does the northwest coast. By 
December the strong gradient regions are similar to those noted in January. 

Overall, the mean of the constructed soil moisture time series seems 
reasonable, with a strong seasonal cycle, although verification is difficult 
due to the lack of extensive soil moisture observations. There are con- 
siderable interannual variations in the soil moisture field, with standard 

deviations around half the magnitude of the soil moisture itself. However, 

* 

this variation is not quite as large as that in the precipitation field, due 
to the cumulative nature of the soil moisture parameter. 

The constructed soil moisture field should be regarded as a first order 
approximation, because a lot of physical details of the soil moisture and 
evapotranspiration calculations have been circumvented. In particular, the 
choice of a single value for the maximum available soil moisture ( 150 mm) 
and the neglect of soil characteristics can be questioned. The use of 
Thornthwaite's empirical formulation for potential evapotranspiration over 
all types of land regardless of vegetation cover is a major simplification 
and is a weak link of the system. We followed this approach because no 
apparent better alternative was available. Nevertheless, the authors believe 
that the constructed fields of soil moisture are suitable as a first order 


approximation 
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Figure Legends 


Fig. 1. 

Fig. 2. 

Fig. 3a. 
Fig. 3&» 
Fig. 4a. 
Fig. 4b. 
Fig. 5a. 
Fig. 5b. 


National Climatic Center (NCC) Local Climatic Data (LCD) station 
sequence numbers. 

Number of complete data years from 1940 to 1980 at each NCC LCD 
station. 

Monthly Mean Station Temperature (Deg. C). 

Standard Deviation of Monthly Temperature (Deg. C). 

Monthly Mean Precipitation (MM/Month). 

Standard Devation of Monthly Precipitation (MM) 

First of the Month Soil Moisture Fraction 

Standard Deviation of the First of the Month Soil Moisture Fraction. 
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Table 1. National Climatic Center Local Climatic Data Stations used. 
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PART II: DEPENDENCE OF LOCAL PRECIPITATION ON ANTECEDENT SOIL 

MOISTURE OVER THE CONTIGUOUS 48 STATES 


ABSTRACT 


In the present study forty years of precipitation data as observed over 
261 Local Climatic Data (LCD) stations in the Continental United States and 
calculated soil moisture fields were utilized to investigate the feedback 
effect of antecedent soil moisture on precipitation. A month-by-month soil 
moisture and evapotranspiration dataset was constructed using observed surface 
temperature and precipitation data in a simple ground wa'ver balance model. 

The monthly precipitation was correlated with antecedent monthly precipitation, 
soil moisture and evapotranspiration separately. The results show a perceptible 
signal indicating a positive feedback between the correlated fields in some 
regions. The maximum positive correlation is found to be in the drought prone 
western Great Plains region during the latter part of summer. There is also 
some negative correlation in coastal regions. The correlations between soil 
moisture, and precipitation, particularly in the latter part of summer, suggest 
that large scale droughts over extended periods may be partially maintained by 
the feedback influence of soil moisture on rainfall. In many other regions the 
lack of positive correlation shows that there is no simple answer such as higher 
land-surface evapotranspiration leads to more precipitation, and points out 
the complexity of the influence of soil moisture on the ensuing precipitation. 
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INTRODUCTION 


Whether changes in land-surface evapotranspiration can significantly 
influence the future precipitation, is a question which has generated some 
controversy in the past. Prior to 1937, it was generally believed that the 
land-surface evapotranspiration had an important Influence on the future 
precipitation. However, beginning with the study by Holzman (1937) several 
investigators (Benton et al ., 1954, Budyko, 1958) presented results contrary 
to this belief. Subsequently, Libby (1959) found that 33% of the rainfall in 
the Mississippi watershed came from local evaporation. For a detailed history 
of this controversy, see Stidd (1968). 

Retrospectively, it can be stated that much of the aforementioned contro- 
versy arose because most of these studies were based on annual averages. 

Since, in the wintertime the rainfall is largely produced by large-scale baro- 
clinic systems and the land-surface evapotranspiration is generally small, 
there is no reason to expect a strong correlation between the two. However, 
the summertime situation is quite different. Generally, over the land surface 
the evapotranspiration exceeds the precipitation, with the exception of 
raonsoonal regions. Also, the precipitation is largely of convective character, 
and derives a significant amount of precipitable water from the mixed layer 
(ML) of the planetary boundary layer (PBL). Therefore, we naturally expect a 
significant correlation between evapotranspiration and precipitation. 

Observationally very little has been done to understand this influence. 

Soil moisture measurements at climatic data stations have been virtually non- 
existent. This has been a great handicap for similar correlation studies 
as well as for general circulation simulation studies. To date, there are only 
a few systematic studies (e.g. Schickendanz , 1976) which shed some light on 
this subject. By employing methods of statistical inference Schickendanz showed 
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that the observed precipitation in the western Great Plains increased after the 
introduction of farming and artificial irrigation. He attributed this finding 
to the influence of increased evapotranspiration. Other researchers have made 
conjectural references to the importance of this influence. For example, 
Sutcliffe (1956) stated that in summertime there is a strong possibility that 
the same moisture can precipitate several times by continual overturning. 

Namias (1959) suggested that the anomalously high rainfall in the summertime 
in Texas in 1957 was perhaps related to the soil moisture anomaly produced by 
an antecedent wet springtime. 

In a recent study of the spatial coherence of monthly precipitation over 
the U. S., Walsh et al. , (1982) found maximum as well as most persistent correla- 
tions between monthly precipitation and the Palmer Drought Index to occur in 
the growing regions of the Midwest. They argued that this was due to the greater 
moisture holding capacity of the soils there. They also noted that the monthly 
precipitation was both least spatially coherent as well as least persistent 
during the summer months. 

In the last ten years, the subject of feedback between land surface 
evapotranspiration and precipitation in numerical models has also drawn consid- 
erable attention. Studies by Walker and Rowntree (1977), Miyakoda et al. 

(1979) and Shukla and Mintz (1982) have indicated the important role of 
soil moisture anomalies on the simulated precipitation. Shukla and Mintz (1982) 
further inferred that the influence of a global scale soil moisture anomaly in 
the summertime can be quite complex. 

In this study, we have examined the relation among monthly precipitation, 
evapotranspiration and soil moisture at 261 Local Climatic Data (LCD) stations 
over the contiguous 48 states in the United States (C48S). The soil moisture 
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and evapotranspiration fields were constructed by a surface water balance 
model using the observed precipitation and surface temperature aB detailed 
in Section X of this report. 


n 


2. PHYSICAL BASIS 

The relationship between monthly evapotranspiratlon and precipitation lor 
a region is simply: 

P *» E + C m (1) 

where P, E and C ra represent monthly mean precipitation, evapotranspiration and 
atmospheric moisture convergence in a region. Obviously, (1) neglects moisture 
storage in the atmosphere. In the summertime, generally speaking, evapotrans- 
piration exceeds precipitation which leads to depletion of soil moisture and 
drier ground by the end of summer. Therefore, if C m were independent of E, it 
would be a simple matter to infer that a reduction in E would naturally lead 
to a corresponding reduction in P. However, there is no justification for 
such an assumption, particularly on the time scale of a month or more. There 
are several other feedbacks on the large scale circulation which play an impor- 

. i 
\ I 

tant role in determining the outcome of these effects. i 

A soil moisture anomaly is directly related to antecedent precipitation j 

anamolies (Walsh et al . , 1982). Because the soil moisture has a time scale of i 

’ ^ 

up to a month or so, even in the summertime, it is plausible to examine these | 

feedback effects on a monthly time scale. Therefore, we propose to use monthly 

. I 

averaged fields for this analysis. Secondly, because the evapotranspiration j 

anomaly caused by a soil moisture anomaly can be expected to influence moisture f 

convergence, it is important to separate these effects. Let us try to explain 
this by an example of near zero evapotranspiration, a drought situation. In i 

such a case precipitation, P, equals the net moisture convergence, C m . Since ■] 

evapotranspiration is a moisture source term, its near absence would lead to ' 

moisture deficiency in the ML. The moist static energy, which essentially I; 

Ij 
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determines the onset of moist-cor 'ection, thereby determining convective 
precipitation is given by: 

h “ CpT + gz + Lq (2) 

Absence of evapotranspiration would reduce the Lq term in the ML. On the other 
hand a corresponding increase in the sensible heat flux must show up in the C p T 
term. The increase in the C p T term will largely, though not entirely, compensate 
the reduction in the Lq term. The difference would presumably be due to the 
higher albedo and enhanced longwave cooling of the dry and hotter ground. There- 
fore, it can be argued that a reduction in evapotranspiration lowers the moist 
static energy of the ML which in tarn could lead to reduced moist convective 
activity. This reduction in evaporation reduces moisture in the ML which then 
leads to raising of the lifting condensation level, reduction in rising motion due 
to reduced condensation heating and lower precipitation, which implies a positive 
feedback. Simultaneously, enhanced moisture convergence produced by a relative 
thermal low induced by sensible heating in a drier region introduces a negative 
feedback. With these two influences operating simultaneously, the overall 
influence of a reduction in evapotranspiration on rainfall is not obvious. 

Nearness of other moisture sources, such as an ocean, may provide moisture to 
the ML for convergence into the dry region. Thus the induced moisture convergence 
may compensate the deficit in moisture produced by lower evaporation. How does 
this local soil moisture anomaly then influence the large scale circulation? 

It is difficult to answer this question intuitively. One approach is to look 
at the observational data to understand how different regions respond to soil 
moisture anomalies. The problem is complex because of the innumerable temporal 
and spatial structures of soil moisture anomalies, as well as the effect of the 
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peculiarities of large scale circulation themselves, having their origin in 
the multiplicity of solutions of the earth-atmosphere system, which in turn 
produce the soil moisture anomalies besides being affected by them. The 
current study attempts to identify regions of positive and negative feedback 
influences so as to understand the role of soil moisture on future precipi- 
tation. However, this study is directly applicable only to summertime in the 
continental United States and the conclusions may not hold for other regions 


or seasons. 



I 



3. SURFACE MOISTURE BALANCE 

Mean monthly Local Climatic Data (LCD) consisting of precipitation and 
surface temperature fields, from 261 stations over the C48S, was obtained from 
the National Climatic Center (NCC) in Asheville, North Carolina, The data 
covers a period of 40 years (1940-1980). Missing data periods were ignored as 
opposed to filling the numbers in by interpolation. Fig. 1 shows the station 
locations with number of years of complete data at each station. Time series 
of potential evapotranspiration, actual evapotranspiration and soil moisture 
were calculated at each station. The details of this calculation are given 


in Section I of this report. We particularly focused on the drought prone 
semi-arid region of the western Great Plains, shown boxed in Fig. 1. This 
subset region as well as the entire C48S were statistically analyzed. 


4, STATISTICAL CALCULATIONS 


The observed precipitation and the constructed fields of antecedent soil 
moisture and evapotranspiration were analyzed to obtain the relationship among 
them* Attention was focused on the entire C48S region and its subset ; (the 
drought prone western Groat Plains, To avoid the feedback of the current 
precipitation upon the calculated evapotranspiration, a psoudo-evnpot cons- 
piration was calculated based on the soil moisture and surface temperature on 
the first of the month. Furthermore, the surface temperature was obtained by 
adding its climatological variation to the antecedent observed value. Thus, 
the caculated antecedent soil moisture and evapotranspiration fields do not 
use the furture precipitation or surface temperature data. This ensures that 
the correlations are truelv between future rainfall and antecedent soil moisture 
or evapotranspiration fields. 

a . Correlation coefficients : 

Three separate temporal correlations were calculated between mean-monthly 
precipitation with a) antecedent mean-monthly precipitation, b) soil moisture at 
the beginning of the month, and (c) monthly evapotranspiration at the beginning 
of the month. Of these only the first correlation is truly between the observed 
fields. We obtained correlation coefficients at each station for each month 
using the 40 year time series. This calculation eliminates non-linear bias 
introduced by spatial averaging and is free of the influences of the seasonal 
cycle. In order to isolate the regions of strong correlations, a two-tailed 
T-test was applied to isolate stations showing a strong signal. The correlation 
of rainfall with antecedent fields (predictors), determined purely from the pot 
history of precipitation, ensures the sanctity of this statistical procedure. 
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b) Contingency tables : 
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This is an alternative approach to correlation analysis which displays the 
structure of the relationship between the analyzed fields. For variables 
showing a skewed (non-normal) distribution this approach is used with better 
success. For each of the two fields to be correlated, the values were divided 
into five categories: very small, small, average, large, and very large. An 

event, e^ is designated to belong to a location i, j (referring to i 11 * 1 
column of j^h row) of this table (i » 1, 5 & j « 1 , 5) depending upon the 
category index of each of the variable!* being correlated. The diagonal of 
events, e-y , from i ® j a 1 to i = j « 5, represents the strength of the 
correlation. In a purely random process the number of events, n^j in a 
location should be calculated from: 

n ij “ N Pi pj (5) 

where N is the total number of events and p^ is the probability of i to be 
associated with any value j and vice versa for p j . To facilitate interpretation 
of the contingency table it is normalized by dividing e^j by njj and then multi- 
plying by 100. Thus a normalized contingency table shows positive (negative) 
correlations for numbers more (less) than 100. 

Obviously, there is no unique way to determine the categorical boundaries 
which divide the rows and columns in a contingency table. The authors tried 
several methods, but found the following to be logical as well as sta i,- Jically 
appropriate. The precipitation at each station for each month was first ordered 
by magnitude, and then divided into categories containing 7,5, 10, 65, 10 and 
7.5 percent of the data, from lowest to highest values respectively. Thus, 
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assuming a normal distribution, these categories correspond to standard devi- 
ation cutoffs at ±2cr and ± lo representing strong, weak or noise level signals. 
Each column in the tables represents a different precipitation category. Thus, 
even though a given table may be constructed by compositing data from many 
stations, the location of an event in the contingency table depends only on 
the variability of the field at that station. For example, entries in the 
driest column truly represent locally dry events. The rows of the tables were 
determined for the variable related to the precipitation, and were filled by 
the same procedure. The contingency tables were tabulated for two regions. 

One region included all 261 stations over the C48S, and the other contained 23 
stations over the drought prone western Great Plains. 

A standard X^-test was conducted to see the strength of the signal in 
these contingency tables. It should be noted that by examining the rows and 
columns constituting the boundaries of a contingency table one can examine the 
correlation of drought (wet season) with the antecedent predictors such as dry 
(wet) soil moisture. 
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5. RESULTS 
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For complete monthly means and standard deviations of temperature, 
precipitation and soil moisture fields over C48S the reader is referred to 
the first section of this report. The mean monthly soil moisture shows a 
distinct seasonal cycle, exhibiting near saturated soil during winter, and 
very dry soil in the latter summer months over most of the stations. The 
frequency distribution of the first of the month soil moisture over C48S, as 
well as over the subset region are shown in Table 1. Over both regions the 
soil moisture is lowest during August, September and October, at which time 
the standard deviation of soil moisture is of the same order as the mean soil 
moisture itself. The mean soil moisture over C48S for the driest month, Sep- 
tember, is shown in Fig. 2a, while the corresponding standard deviations are 
shown ir. Fig. 2b. 

The correlation coefficients between precipitation and the generated hydro- 
logy fields over the C48S follow a pattern very similar to that of precipitation 
with itself with a one month lag. Over most of the C48S there is no coherent 
pattern of significant correlations between these fields. September and October 
show the largest number of significant correlations, and have a coherent pattern 
of significant correlations over the western Great Plains. Fig. 3a shows the 
September - October precipitation correlations, while Figs. 3b and 3c show the 
soil moisture - precipitation and evapotranspiration - precipitation correlations 
for October. In these figures correlation coefficients significant at the 80% 
level or above are enlarged, while those at a few stations with less than 8 
years of data are omitted. Also notable in these figures is the coherent 
region of negative correlations along the northwest and the Gulf coasts. 

The September correlations (not shown) are similar, with the western Great 
Plains again being the only coherent region of significant correlations. 
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The contingency tables over CA8S for interdependence of precipitation with 
itself at a one month lag* and for soil moisture and evapotranspiration with 
precipitation are shown in Tables 2a, b and c. The contingency table at the 
bottom right is a composite of July, August, September and October, Ghi-aguared 
test values are included at the bottom of each contingency table. The cor- 
responding tables for the subset region are shown in Tables 3a, b and c. A 
diagonal dominance from upper left to lower right can be observed in the latter 
summer months, particularly September and October, in all six sets of tables. 

This dominance, which represents well correlated fields is stronger over the 
subset region. The July through October composite contingency tables display 
this dominance well, particularly in the two corners, which represent well cor- 
related extreme cases. The large X 2 values reflect the strength of the signal, 
figure A shows the normalised table values in the driest (solid) and wettest 
(dashed) soil moisture categories for the July to October composite ease. Over 
both regions (Pig. Aa and Ab) the same trends appear, although they are stronger 
over the subset region. Dry (wet) soil corresponds to less (more) precipitation. 
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We have isolated some regions and a few months of a significant positive 
relationship between precipitation and antecedent soil moisture. The correla- 
tions over the drought-prone western Great Plains were perceptibly stronger 
than those over the rest of the United States. In coastal regions, particularly 
the west coast, the correlations, though small* often had a coherent negative 
pattern. It should be pointed out that smaller negative correlation coefficients 
are to be expected in negative feedback situations, and their significance 
must be interpreted in light of this. Since recent simulation studies by 
Shukla and Mints (1982) and Sud and Fennessy (1983) with the GLAS GCM have 
also shown that soil moisture influences on precipitation can be both positive 
or negative , the present results are not surprising. Increased moisture con- 
vergence due to higher ground temperature could explain this negative feedback 
in regions near moisture sources, e.g. the west coast of United States. The 
correlation coefficients generally are not y.*.-> large, but nonetheless are 
significant in some regions . The overall distribution of correlation coeffi- 
cients is somewhat complicated, which indicates the complexity of the feedback, 
as well as the strong influence of many other factors not taken into account. 
However, the summertime contingency tables do show that for extreme cases, 
particularly for very dry events , there is a clear dependence of precipitation 
on soil moisture. The July to August correlations and the August contingency 
tables seem to be an exception to this rule, showing no clear-cut relation 
between these variables. This could perhaps be due to the effect of local 
irrigation, which Schickedanz (1976) has shown to have a significant effect on 
the local precipitation. Schickedanz found increased local precipitation in 
irrigated regions of the western Great Plains during June , July and August, 
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with the strongest effect occurring during July. As we do not account for 
thi6 effect in the current study, it would lead to a degradation of our results. 

The correlations were obtained by constructing soil moisture and evapo- 
transpiration fields from a simple water balance model, which itself may only 
be adequate in semi-arid regions. Neither the field capacity (w* 5 150 mm ) 
nor the use of Nappo's (1975) evapotranspiratlon function are so suitable in 
forested regions found over many parts of the United States. We must also 
point out yet another drawback of a study such as this. The actual observations 
have many other influences such as changes in the boundary forcing and Initial 
conditions in the real earth-atmosphere system. Some of the competing influences, 
6uch as sea-surface temperature anomalies, the initial state of the atmosphere, 
etc. are known to be important, but are not a part of this analysis. Perhaps, 
the true magnitude of the impact of soil moisture anomalies cm the precipitation 
can only be studied by controlled experiments with GCMs. 

For the present study, it is important to point out that some statistically 
significant Influence has been found, both in space and time, between precipi- 
tation and surface evapotranspiration parameters. This influence may perhaps 
play a key role in determining the late summer precipitation in such regions and 
may therefore have an important role in the maintenance of droughts. 
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Fig. i. Number of complete data years from 1940 to 1980 for NCC LCD U.S. 
stations. 

Fig. 2a. Mean September soil moisture fraction over NCC LCD U.S, stations. 

Fig. 2b. Standard deviation of September soil moisture fraction over NCC LCD 
U.S. stations. 

Fig. 3a. September to October precipitation correlation coefficients. 

Fig. 3b. October first of the month soil moisture to October precipitation 
correlation coefficients. 

Fig. 3c. October first of the month evapotranspiration to October precipitation 
correlation coefficients. 

Fig. 4a. Graphical display of the normalized probability of correlation between 
antecedent soil moisture and precipitation over C48S. 

Fig. 4b. Graphical display of the normalized probability of correlation between 
antecedent soil moisture and precipitation over subset region. 
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TABLE 1 

SOIL MOISTURE FREQUENCY DISTRIBUTION 
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TABLE IX. A W V 

NORMALIZED CONTINGENCY TABLES FOR C48S 
MONTHLY PRECIPITATION WITH ANTECEDENT PRECIPITATION 
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NORMALIZE!) CONTINGENCY TABLES FOR C48S 
MONTHLY PRECIPITATION WITH ANTECEDENT SOIL MOISTURE 
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NORMALIZED CONTINGENCY TABLES FOR SUBSET REGION 
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MONTHLY CONTINGENCY TABLES FOR SUBSET REGION 
monthly precipitation with ANTECEDENT EVAPOTRANSPIRATION 
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